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We discuss the sensitivity of the three-particle azimuthal cumulant method for a search and study 
of conical emission in central relativistic A + A collisions. Our study is based on a multicomponent 
Monte Carlo model which include flow background, Gaussian mono-jets, jet-flow, and Gaussian 
conical signals. We find the observation of conical emission is hindered by the presence of flow har- 
monics of fourth order ( 114) but remains feasible even in the presence of a substantial background. 
We consider the use of probability cumulants for the suppression of 2 nd order flow harmonics. We 
' find that while probability cumulant significantly reduce v\ contributions, they also complicate the 

cumulant of jets, and conical emission. The use of probability cumulants is therefore not particu- 
f"^ , larly advantageous in searches for conical emission. We find the sensitivity of the (density) cumulant 

£SJ ' method depends inextricably on strengths of V2, V4, background and non-Poisson character of par- 

ticle production. It thus cannot be expressed in a simple form, and without specific assumptions 
about the values of these parameters. 
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0) ; I. INTRODUCTION 



Observations of away-side dip structures in two-particle correlations measured in Au + Awcollisions at yjs nn = 200 
GeV have stimulated renewed interest in the notion of conical emission. The passage of a parton through dense matter 
at speeds greater than the speed of sound is predicted to lead to the production of a Mach shock wake resulting in 
conical emission pattern that may explain the observed two-particle correlations [j], El, El, 0, IE IE 0, IE 13, E3, EH- 
Identification of Mach cone is of great interest because it could provide an experimental determination of the speed 
of sound in the dense medium produced in high energy A + A collisions [3, [l2| . A sonic boom is also expected for a 
heavy quark propagation based on Ads/CFT calculation [l3|. Cerenkov radiation produced by a superluminal parton 
. traversing a dense medium is expected to generate a similar signature fLU [TBL [Tsl. fl7l [l8j . However, other production 
mechanisms have been proposed to explain the two-particle correlation data. These include large angle gluon radiation 
0,113], Path length effects 0], collective flow, and jet deflection EE EI EE EE ■ 

Three-particle correlation measurements were proposed to gain further insight in the particle production mechanism 
00 ! leading to the away-side deep structure seen in two-particle correlations. Various methods have been suggested 
■ and are currently pursued to carry such analyses. A number of these analyses are carried out using a flow + 
signal decomposition, i.e. an adhoc flow background is assumed and subtracted based on the ZYAM approximation 
H, EE E3, EH, EE EE EE EE EE EE EH- Note however that, as Borghini pointed out, momentum conservation 
can have a significant impact on two- and three-particle correlations [3^ ] . The magnitude of momentum conservation 
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effects is however difficult to estimate as discussed recently by Chajecki et al. [40[ ■ A cumulant based method was 
proposed by the author 0] and is also used by the STAR collaboration 0] to search for conical emission. This 
method presents the advantage that the extraction of a three-particle signal does not require any model assumptions. 
Obviously, one can also carry a model based decomposition of the extracted three-particle signal, i.e. the cumulant, 
to estimate its various components, and seek evidence for conical emission. 



In this work, we further study the properties of the three-particle cumulant observable first proposed in Ref. 41] 
for a search for conical emission, and studies of particle production dynamics. Three particles correlations, designed 
to search for conical emission, use a high p t particle as a jet tag. The azimuthal angle of this high p t particle should 
approximately correspond to the direction of the parton initiating the tagged jet. It is further assumed the emission 
of this parton is surface biased, and directed outward, approximately normal to the surface of the medium. The 
direction of the high-pt particle thus provides a reference to study the propagation of the away-side parton initiating 
the second jet. Estimates of the sound velocity in the quark gluon plasma suggest the wake produced by a high energy 
arton should produce particle emission at a Mach angle of the order of one radian relative to the parton direction 
This corresponds to enhanced particle emission at angles of order 2 radian with respect to the high-p f jet tag. An 
unsuppressed away-side jet should on the other hand lead to particle emission at an angle of approximately tt relative 
to the jet tag. Deflection or radial flow effects should produce a broadening of the away-side jet relative to the high-p t 
tag. 
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We define the notation and variables used in this paper in Section|ITJ We next discuss, in Section ITTTI simple models 
of particle production including di-jets, flow, conical emission, as well as the effects of differential jet quenching, 
hereafter called jet-flow. We show in Sec. IIIIBI that in the case of Poisson statistical particle production, the three 
particle cumulant (hereafter noted 3-cumulant) associated with azimuthal anisotropy reduces to a simple expression 
involving only non-diagonal Fourier terms, fn view of the large values of elliptic flow v 2 observed in A + A collisions 
at RHIC, one expects the leading non-diagonal term should be of order V2V2Vi. However, one finds experimentally 
that fluctuations in the number of produced particles are non-Poissonian. We thus consider in Sec. IIVI the impact 
of such non-Poissonian fluctuations, and find they imply the presence of sizable v\ terms in the threc-cumulant that 
may complicate the observation of conical emission signals. We show the v\ terms may however be suppressed if one 
adopts a modified version of the cumulant based on probability densities rather than number densities, and suggest 
this modified cumulant as an alternative means of studying three-particle correlations. 

The cumulant method is a relatively complicated analysis technique, which may in principle be sensitive to various 
instrumental effects. Wc discuss a specific implementation of the method based on normalization to single particle 
distributions in Sec. [V] We show, based on the conical emission model introduced in Sec. I VII and an adhoc 
parameterization of the detection efficiency, that normalization by single particle distributions leads to a robust 
analysis technique. 

A key issue in the search for conical emission is the sensitivity of the method used to extract the three- particle 
correlation information. We discuss, in Section [VTl the sensitivity of the cumulant method based on a simple jet and 
conical emission Monte-Carlo model. 

Our conclusions are summarized in Section [VIII 



II. DEFINITIONS AND NOTATION 

The conical emission search method introduced in Ref. [4l| is based on the observation of three-particle densities as 
a function of the relative azimuthal angles between the measured particles. The presence of three-particle correlations, 
and possible signal for conical emission is extracted using the 3-cumulants. A high p t particle is used as a jet tag, 
and proxy for determination of the direction of the jet. Two lower p t particles are used to probe the structure of the 
near-side jet (i.e. the jet singled out by the high-p t tag particle) and search for conical emission on the away-side. The 
jet tag is herein referred as particle 1 while the two associates are labeled as particles 2 and 3. The three particles 
are detected in the collision transverse plane at angles tpi, tf2, and 933 relative to some arbitrary reference frame. The 
single, two-, and three- particle densities are noted as follows: 

pi(tp%) = dNi/d(pi (1) 
P2(<Pi,<fj) = dN 2 /dipidtpj 
P3(<Pi,<Pj,fk) = dN 3 /difidipjdipk 

where the indices i, j, and k take values 1, 2, and 3. 

The three-particle or triplet density, ps, corresponds to the average number of particle triplets observed per collision. 
The particles of a given triplet are however not necessarily correlated. Indeed, the three particles may originate from 
one, two, or three (distinct or not) production processes (e.g. radial flow, elliptic flow, resonance decay, jets, Mach 
cone, etc) . Cumulants are designed to extract the three-particle correlation c omp onent from the three-particle density. 
They were first discussed in the context of particle physics (see for instance [43[ 53 ) and are now used in a variety of 
analyses. It is the purpose of this work to study specific aspects of the cumulant method introduced in [4l| for studies 
of the shape and strength of the signal expected from different types of processes and to characterize the robustness 
and sensitivity of the method. 

The 2- and 3-cumulants C 2 and C3 arc defined as follows: 

C 2 (ipi,ipj) = p2(<Pi,tPj) ~ Pi(<Pi)pi(<Pj) 
C 3 (<Pi,<Pj,<Pk) = P3(<Pi,<Pj,<Pk) ~ P2(<Pi,<Pj)pi(tpk) (2) 
~ P2(<Ph<Pk)Pl(<Pj) ~ P2(<Pj,<fik)Pl(<Pi) 

+ 2pi{<Pi)Pi(<Pj)pi(<Pk) 

As described in Ref. [4l[ the two-, three-particle densities and cumulants are straightforwardly corrected for detector 
inefficiencies provided the three- and two-particle detection efficiencies may be factorized as products of respectively 
two and three single particle efficiencies. The robustness of this correction procedure is discussed in Sec. |V] on the 
basis of Monte Carlo simulations. 
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Correlation functions in terms of relative angles, C 3 (Atpij , A(pi k ) , are formally obtained by integration of the 
cumulants C 3 (tpi,tpj, ipk) with constraints Aifj,j = tp^ — ipj. 

C 3 (Atpij , Atpik) = I ' C 3 (<pi,<pj,<pk) (3) 

x 5(Aipij - <fi + ifj) 

x 8(Aip tk - tpi + ip k )dipid<pjd(p k 

In practice, this is accomplished by binning C 3 ((fi, ipj, ip k ) and C 3 (Atpij, Aip ik ) into arrays (e.g. 72x72x72 and 72x72 
respectively) and summing the elements of C 3 (tfi, ifj, ip k ) to obtain C 3 (A(fij, A(fik) as follows. 

72 

C 3 (p,«)= E C 3 (i,j,k)S(p-i + k)6(q-i + k) (4) 

i,j,k— 1 

for p=l, ■■-,72; 5=1, ...,72. Given the finite statistics, and the large memory requirement implied by three-dimensional 
arrays used in this type of analysis, care must be taken when binning the densities. We found the use of 72 bins, 
for analysis of data samples collected e.g. by the STAR experiment [45|, HI] enables sufficient angular resolution and 
statistical accuracy. 



III. SIGNAL MODELING 



We use specific models for jet, di-jet production, collective flow, and conical emission signals in order to assess 
the sensitivity of the cumulant method for a search of conical emission. The models are described in the following 
sub-sections while our study of the sensitivity of the cumulant method is presented in Section IVII 



A. Di-Jet Production 



Jet production is characterized by the emission of particles in a cone (in momentum space) centered on the direction 
of the parton that produces the jet. We consider di-jet production restricted to central rapidities (near 90° relative 
to the beam direction), and assume the number of di-jet per event in the acceptance of the detector, J, varies event- 
by-event. The intra-jet particle multiplicity depends on the jet energy but is order three (3) at RHIC energy [47j . 
We consider back-to-back jets in azimuth, but possibly different longitudinal momenta. We denote the number of 
jet particles measured within each kinematics cut "i" as Aj. and A! i respectively for the near and away side jets. 
The emission of particles relative to the jet axis is described by a probability distribution P aS soc(Q), where 6 is the 
angle between the particle momentum and the momentum vector of the parton originating the jet. Projection of this 
distribution in the transverse plane leads to a probability distribution, Pj(ip), function of the azimuthal angle, <p, 
between the parton and the particle direction. Following (4l| . we further simplify the jet model and use a Gaussian 

azimuthal profile, G 2 ((fii — <fi a ;<Ti) with G 2 (x\a) = (y/2na) exp (— x 2 /2a 2 ) where <fi, cf> a , and a% are the emission 
angle (in the lab frame) of measured particles, the emission angle of the parton, and the width of the jet, respectively 
For this illustrative model, we assume the emission angle and the jet multiplicity are not correlated to the collision 
reaction plane. We further assume one can decouple the number of particles, in the measured kinematic range of 
interest, from the jet profile. The number of jets, their multiplicity, and directions are not known, and must therefore 
be averaged out. The 3-particle jet cumulant is given by the following expression. 



C 3 Jet (A^ 12 ,A^ 13 ) = (2tt)- 1 (J) 



(A X A 2 A 3 ) G 3 (Aip 12 , Aip 13 ; a 1 ,a 2 ,cj 3 ) 
+ (A^A's) G 3 (A Vl2 ,Aip 13 - tt; a u <r 2 , a' 3 ) 
+ {AxA' 2 A 3 ) G 3 (Aip 12 - tt, A^ 13 ; a lt a' 2 , a 3 ) 
+ (AiA^A'z) G 3 (Aip 12 - tt, A^ 13 - tt; ax, a' 2 , a' 3 ) 

{A X A 2 



(5) 



+ (2n)- 2 ((J(J -!))-< J) 2 



(A 3 +A' 3 ) G 2 (A<p 12 ;a 12 ) 



+ (A 2 A 3 ) 
+ (AiA> 2 ) 
+ (AiA> 3 ) 
+ (A' 2 A> 3 ) 



(A 2 + A' 2 }G 2 (A(p 13 ;a 13 ) 
(Ax) G 2 (Aip 13 - Aip 12 ;a 23 ) 



(As + A' 3 )G 2 (A<p 12 - tt; 
+ A' 2 ) G 2 (Aip 13 - n; 
G 2 (Aipi 3 - Aip 12 - 



(A 2 
(Ax) 



'13) 



J 23) 



+ (27T)- 3 (<J(J - 1)(J - 2)) - 3 (J(J - 1)) (J) + 2 (J) 3 ) (Ax) (A 2 + A' 2 ) (A 3 + A' 3 ] 
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with 

G3(x 1 ,X2;a 1 ,a2,a- 3 ) = {2ny 1 a^ 3 X exp(-(<7fxi + 0-2X2 + 0i(a:i - x 2 ) 2 ) / 2af <23 ) (6) 

where Xij = Xi — Xj, afj = erf + cx| and crfj k = <J 2 o 2 + of cr^ + o 2 o\- Note G3 is defined such its integral over x\ and 
X2 is unity. The coefficient (J) determines the average number of jets per event. Coefficients A 2 , A3 describe the 
number of particles associated with the near side jet. while A[, A 2 , A' 3 correspond to the number of particles associated 
with the away side jet observed within the kinematic cuts used for the study of the correlation functions. We note 
that if the number of jets. J, in each event is determined by a Poisson process, then one has (J(J- 1)) = (J) and 
(J(J — 1)(J — 2)) — 3 (J (J — 1)) — 2 (J) = 0. The constant term, and the terms containing a 2-particle dependence 
in G 2 then vanish in the above expression. The jet 3-cumulant is plotted in Figure 1, in arbitrary amplitude units, 
for (7i = 0.15, <7 2 = 03 = 0.2, a' 2 = <j' 3 = 0.4, (J) = 1, Ai = 1, A 2 = A3 = 2, and A' 2 = A' 3 = 1. We further assume 
for simplicity that (AiAj) = (Ai) (Aj) and (A\A 2 A3) = (A\) (A 2 ) (A3) although it is unlikely realized in practice. 
Panel (a) presents the cumulant for the Poissonian case, whereas panel (b) displays a non-Poissonian case for which 

( J( J — 1)) — (J) 2 J (J) 2 = 0.05. The peak centered at Aipi 2 = Aipi3 = corresponds to two particles associated 

with the high p t particle, and is referred to as the near side jet peak. Secondary peaks at (0, 7r), (tt, 0), and (jt,tt) 
correspond to one or both the associates being detected on the away side. The bands seen in Figure Q] (b) stem from 
the non-zero two-particle contributions to the 3-cumulant for non-Poissonian events. Various effects may alter the 
strength and shape of the jet correlations. Interactions of the parton with the medium may produce jet broadening, 
and deflection. Given jet emission is expected to be surface biased in A + A collisions (in view of recent measurements 
of small RaA: and the disappearance of the away side jet (48j). one expects the near side peak to have similar width in 
A + A collisions as in p +p collisions, while the away-side peak should be broader. One might also observe additional 
broadening along A(p± 2 , and Aipi3 due to parton scattering, and radial flow. 




FIG. 1: 3-Cumulants calculated based on the Gaussian jet model discussed in the text, (a) Poisson case, (b) Non-Poisson case 
with <J(J - 1)) - (J) 2 = 0.05(J) 



B. Collective Flow 

Flow, or collective motion, is an important feature of heavy ion collisions at relativistic energies. It manifests itself 
by a modification of transverse momentum (p t ) spectra relative to those observed in p + p collisions, and by azimuthal 
anisotropy of produced particles. In this section, we focus on azimuthal anisotropy arising in non-central heavy ion 
collisions. We decompose the azimuthal anisotropy in terms of harmonics relative to an assumed reaction plane. The 
conditional probability, Pp(ipi\ip), to observe a particle at a given azimuthal angle, tpi, with the reaction plane at 
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angle, ip, is written as a Fourier series. 



(7) 



where 



F 2 (A(p;Vr, 



= 1 



cos (m (Aiyj)) 



(8) 



The Fourier coefficients v m (i) measure the m-th order anisotropy for particles emitted in a selected kinematic range i. 
Measurements have shown the second order (elliptical) anisotropy can be rather large in Au + Au collisions at RHIC 
while first and fourth order harmonics are typically much smaller. STAR measurements show the fourth harmonic 
scales roughly as the square of the second order harmonic (vi ~ l.luf) [48|,[4!|. The third and fifth harmonics are by 
symmetry null at r\ = for a colliding system such as Au + Au, and expected to be rather small at other rapidities. 
Higher harmonics are most likely negligible. We neglect possible event-to-event fluctuations of these coefficients given 
it emerges from recent works disentangling flow fluctuations and non-flow correlations is difficult. We describe the 
probability of finding N particles in the kinematical range i according to probability P{N). The exact form of this 
probability is not required. Only its first, second, and third moments are needed. The joint probability of measuring 
Ni particles at an angle <pi while the reaction plane angle is at ip is given by Pp{(fi,Ni, ip) = PF{(pi\tp)PF{Ni)Pnp(ip) 
where Prp{iP) = (2tt)~ 1 is the probability of finding the reaction plane at a given angle ip. Integration yields the 
single particle density pi(fi) = (27r) _1 {Ni). The flow 2- and 3-cumulants are given respectively by (see Ref.pH for 
a derivation of these expressions): 



C 2 F (A^y) = (2n)- 2 {{N i N j )F 2 (A< Pij ;v m {i)v m {j))-{N i ){N j )) 
= {2n)- 2 {N i N j ) 

x ^1 - dij + 2 v m {i)Vm(J) cos (mAy)jj)j 



(9) 



and 



Of (</>,, <fij, tp k ) = (2tt)- 3 (NiN 2 N s ) [M<Pi,<Pj,<Pk) 

+ (l-/y fc )* 2 (Apy) 

+ (l-/jjy)*2(Apik) 

+ (l-/ifci)*2(A^- fc ) 

H~ 1 fijk fikj fjki ^Qijk] 



(10) 



where 



(11) 



m— 1 



and 



$3(</>i> Vji^fc) = 2 ^2 vi iP v 2>m v 3 , n 



J p,m-\-n 



cos {pipi — nupj — nipk) 



+i5m, P +n cos (— ptfi + rrupj - nipk) 
+S n ,m+ P cos {-pipi - rnipj + nipk) 



(12) 



Vi }m correspond to the anisotropy coefficients of order m for particle i. The coefficients dij, fijk, and gijk are defined 
as follows 



_ {N) {Nj) 
ij {NiNj) 



fijk 
9ij k 



{NjNj) (N k ) 
{NiNjNk) 
{Nj) (Nj) {N k ) 
{NNjN k ) 



(13) 
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FIG. 2: 3-cumulant flow signal calculated using (a) only V2V2V4 terms, (b) V2V2V4 and V2V2+V4V4 terms in 1:0.5 ratio, (c) V2V2V4 
and V2V2 + V4V4 terms in a 1:2 ratio, (d) only V2V2 + V4V4 terms. The V2V2V4 and V2V2 + V4V4 terms are calculated according to 
Eqs. I1H14I See text for details. 



We find that the flow 3-cumulant involves, in general, an arbitrary combination of second (e.g. uf) an d third order 
(e.g. V2V2V4) terms. If particle production is strictly Poissonian, the coefficients fijk and gijk are unity. The second 
order terms of the cumulant thus vanish, and the 3-cumulant (flow only) then reduces to: 



CF,Poisson ( \ 
3 {<Pl,V2,¥>3) 



(2tt) 3 <JViJV 2 JV 3 ) #3(^1,^2,^3; w P (l)«m(2K(3)) 



(14) 



which features only non-diagonal terms v m {i)v n {j)v p {k) 1 i.e. with m ^= n, m ^ p, or n 7^ p. This is illustrated in 
Figure [2] (a) which shows the flow 3-cumulant calculated using only V2V2V4 terms. In general, particle production is 
non-Poissonian, and second order terms must be considered. The flow 3-cumulant shown in Fig. [2] (b) and (c) includes 
v\ and v\ added in ratios 1:0.5, and 1:2 respectively to the V2V2V4 terms. Only v\ and v\ terms are included in Fig 
[5] (d) illustrating a case where non-Poisson fluctuations are very large, and thereby dominate the 3-cumulant. The 
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shape of the 3-cumulant thus depends significantly on the strength of the non-Poisson (2 nd order terms), as well as the 
magnitude of the flow coefficients. The interpretation of measured 3-cumulants may thus be considerably complicated 
by the presence of non-Poisson fluctuations. 

Experimentally, one can estimate the magnitude of second order (diagonal) terms from measured total numbers of 
triplet, pair, and single particle densities. STAR observes based on data from Au + Au collisions at 200 GeV [46|, the 
coefficients /123 and /231 deviate from unity. The size of the deviation scales qualitatively as the inverse of the event 
multiplicity. STAR also finds the magnitude of the deviation depends on the specific particle kinematic ranges used to 
carry the analysis. The particle production process is manifestly non-Poissonian. Since the bulk of produced particles 
exhibits elliptic flow, this implies the conditions fijk = gijt = 1 are not verified in practice. The second order and 
constant terms of the cumulant do not vanish, and may in fact be sizable. The 2 nd second order terms result from 
number fluctuations, and may in principle be estimated (and thus fitted) on the basis of measurements of fijk- Note 
however that multiplicity fluctuations occur for flow-like processes as well as all other types of production processes 
(including jets, conical emission, etc). It is thus non-trivial to unambiguously evaluate the proper magnitude of the 
diagonal flow terms based on ratios of measured yields. 

Effects of flow fluctuations are not specifically addressed in the above flow model. Note however that, as defined, 
the 3-cumulant measures averages of (v^), and (viVjVk) coefficients rather than averages of (v n ). The 3-cumulant 
thus implicitcly features fluctuations and non-flow effects. Given it is difficult to experimentally distinguish flow 
fluctuations and non-flow effects (see for instance the review by Voloshin et al. (50|). analyses attempting explicit 
subtraction of flow contributions (from three particle densities) based on the above equations, are thus intrinsically 
non-robust. 



C. Conical Emission 

Mach cone emission of particles by partons propagating through dense QGP matter was proposed by Stoecker Q 
to explain the peculiar dip structure found at 180° in two-particle correlations reported by the STAR and PHENIX 
collaborations, and is the subject of many recent theoretical investigations [l], 0, y, 0, H, 0, 0, H, S 0, EH, HH HE HH 
I2H [27j . The concept of Mach cone emission is based on the notion that high momentum partons propagating through 
a dense QGP interact with the medium and loose energy (and momentum) at a finite rate. The release of energy 
engenders a wake that propagates at a characteristic angle, the Mach angle, determined by the sound velocity in the 
medium. Authors of Ref. (jj estimated the speed of sound in the QGP to be of the order of v s w cf HIC w 0.33. The 
Mach angle should thus be of the order of 70° relative to the away-side parton direction. We use this prediction to 
motivate a simple geometrical model of conical emission. The near- side jet is described using a Gaussian azimuthal 
profile. Away side particles are emitted at a cone angle of 1 radian with a scatter of 0.05 radian relative to the 
away-side direction. Figures [7] (a-b) present the 3-cumulant obtained with a Monte Carlo simulation based on two 
million events. Given that Mach cone particles are produced at 1 radian from the away-side direction and roughly 
normal to the beam direction, narrow Jacobian peaks are seen in the three-particle correlations. A strong dip is 
present at 180° in the two-particle correlations, while in the three-particle correlation a clear spacing is found between 
the peaks. In this simple model, the finite width of the peaks is due in part to the finite width of the trigger jet, and 
in part to the scatter imparted to the cone particles relative to the away-side direction. In practice, one might expect 
additional broadening of the cone because the speed of sound changes through the life of the QGP medium, and given 
the finite size of the medium. Radial flow might also significantly alter the correlation functions. As discussed in 
[4l| , the details of the correlation shapes clearly depend on assumptions made about the kinematics of the away-side 
parton. 

D. Flow x Jet Correlations 

The measured nuclear modification factor, Raa, defined as the ratio single particle yields measured in A + A 
collision to those measured in p + p interactions (scaled by the number of binary collisions in A + A) sug gest the 
propagation of jets through the medium is severely quenched, and therefore features a surface emission bias [48j. In 
this scenario, partons propagating through the medium loose a large fraction of their energy. This results into jet lost 
or lower energy jets with smaller particle multiplicity. The path length through the medium determines the amount 
of energy loss and quenching. At a given point of emission, near the medium surface, a parton propagating outward 
in a direction normal to the surface should have shortest pathway through the medium, and minimal energy loss. 
Partons produced at the same location, but emitted in other directions would have longer path through the medium 
and suffer larger energy loss. In this scenario, one thus expects jet yields (and consequently the high p t trigger) to 
be correlated to the reaction plane orientation. Neglecting disturbances imparted to the medium by the propagation 
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of the jet (or parton), wc model the jet dependency on azimuthal angle relative to the reaction plane with Fourier 
decomposition. Specifically, we write the probability of the jet being emitted at angle (f> while the reaction plane is at 
angle ip as 

P(0|V>) = l + 2^a„cos(n(0- ^)) (15) 

n 

where the coefficients a n represent the effect of the differential azimuthal attenuation. We parameterize the jet 
multiplicity and azimuthal width using associated yields, and Gaussian widths that do not depend on the azimuthal 
direction. As in Sec. IIII A( the average number of jets per event is written (J) while the number of particles 
associated with a given jet is (A)- We also assume the presence of a flowing background consisting on average of 
(Bi) particles. The single particle yield is thus pi,j®F{<fi) = (27r) _1 ((J) (A;) + (Bi)). Given the production of both 
jets, and background particles is correlated to the reaction plane, one ends up with flow-induced correlations between 
all particles, even those produced by the jets. The two-particle cumulant involves three sets of terms, namely jet-jet 
(J® J), jet-background (J®B), and background-background (B®B). Similarly, the three-particle cumulant involves 
J ® J ® J, J ® J ®B, J® B® B, and B® B® B terms. 

A jet-flow model was already discussed in [4lj . We here re-derive, and express the various terms of the model with 
a more compact and intuitive notation. The terms, B ® B and B ® B ® B are identical in form to those considered 
in Sect IIII B[ We need to discuss only the terms J ® J, J ® B, J ® J ® J, J ® J ® B, and J ® B ® B. We begin 
with J ® J, and J ® J ® J. These terms contain contributions where two (three) of the particles are part of the 
same jet, and others where one, or two particles are not from the same jet. When particles are from the same jet, the 
correlation to the reaction plane is unimportant (to the extent, in our model, that the width of the jets do not change 
with their orientation relative to the reaction plane). The corresponding terms are thus identical to those obtained in 
the absence of flow. The two-particle density includes a term with two particles from the same jet (identical to that 
discussed in Sec. IIII A|) . and a term with particles from two different jets. This last term is given by: 



(J(J - 1)) (A) (Aj) J dil)d4 a G{ipi - <t> a ) (16) 
xP(0 a |V) /<%G(^ - ^)P(^|V) 



Integration, and inclusion of the first term yields the jet component of the two-particle density. 

P2 et ° nly (<Pi,<Pj) = (J) (AiAj) G 2 ((fi - fj;o- 12 ) (17) 
+ (J (J - 1)) (A) (Aj) F 2 (cpi - W <»n(i, j)) 

where o\ 2 = a\ + erf, CL n (i,j) = v' n (i)v' n (j), with v' n (i) = a n exp(— n 2 of/2). Similarly, the three-particle density 
includes a term with three particle from the same jet, two particles from the same jet, and another with all three 
particles from different jets. The J ® J, and J ® J ® J part og the cumulant is as follows. 

p( etonv (ipi,(pj,(p k ) = (J) (Aj A,Afc) G$(ipi, ifj, ip k ; Oi, aj, a k ) (18) 

t a 2 -2 \ 



+ (J(J - 1)) (A A,-) (A k ) G 2 (A^)F 2 + -^-ifj - <p k ;i%{k) + (jki) 



+ (J(J - 1)(J - 2)) (A) (Aj) (A k ) N 3 {<pi, cpj, Vk;<(i)y m (j)y p {k)) 

The functions G3, F3, F 2 were defined in previous sections. The symbol (jki) indicates identical terms obtained by 
permutation of the indices. The term G3 is identical to that obtained for non-flowing jets. The term of the last row 
corresponds to a flow-like correlation between the three measured particles, albeit with a strength which depends 
on the fluctuations of the number of jets per event, and the jet associated multiplicities in the kinematic ranges 
considered. The second row contains three similar terms that embody the jet-flow correlations: two particles are 
correlated because they belong to the same jet, while the third is correlated to the first two because of reaction plane 
dependencies. One notes that in the above expression, F 2 depends on the angle of emission of all three particles (at 
variance with assumptions made in Ref. [3lL l32l] ) in a non trivial way. Indeed, in general, the jet correlation widths o~i 
depend on the particle momentum ranges: the ratios of /afj are arbitrary (non-integer) values. F 2 therefore involves 
inharmonic flow components. This implies it is inappropriate to model jet-flow cross-terms as a simple product of 
jet-like correlations, and flow- like correlations, as in Ref. [U, H^, for purposes of background subtraction. 

We next consider the component J ® J ® B, and J ® B ® B of the 3-cumulant. Following the steps used for the 
derivation of the pure jet components, one finds it also includes a term proportional to F 2 G 2 indistinguishable from 
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J ® J Cg> J. One has 



(J) (AiA 2 ) (B 3 ) G 2 (A<p 12 ; a 12 )F -§- Vl + -±<p 2 - <p 3 ; <(3) (19) 

\ tT 12 a 12 / 

The J <g) B (g) B compoment of the 3-cumulant contains a non-diagonal flow term as follows 

(J) (Ax) (B 2 B 3 )F 3 (ip i , 9j ,ip k ;v' n (i),v' m (j),v' p (k)) (20) 

This term has the same functional dependence and is indisguishable from the term on the third line of Eq 1191 
Assemblinging all components, one obtains the JET <g> FLOW 3-cumulant: 



C 3® F ( ( Pu ( Pj> ( Pk) = (J) {AiAjA k )G 3 {}Pi,<pj,ip k ;<Ti,(jj,o k ) 

+ ((J(J - 1)> - (J) 2 ) (AiAj) (A k ) G 2 (Acp t3 ) + (jki) 



(J 2 G 2 \ 

+ (AiAj) ((J(J - 1)) (A k ) + (J) (B k )) G(A i p ij ;a ij )^ 2 I -%-<Pi - -j-¥>i - fk'X(k) + (jki^l) 
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- (J (J - 1)) (Aj) (Aj) (A k ) F 2 (t Pi - tpjioniij)) + (jki) 

+ (J(J - 1)(J - 2)) (At) (A,) (A k ) AT 3 fe, <p k ; v' n (i), v' m (j), v' p (k) 

+ (J) (Ai) (BjB^N^ipi^^ipk-X^Mjlv 1 ^ 

- (B t Bj) (J) (A k ) F 2 fa - <pj;v n (i)v n (j)) + (jki) 
+ (2ir)~ 3 (BiBjB k ) $ 3 (<pi, ifj , ip k ;v p (i)v m (j)v n (k) k ) 

+ (2iry 3 (B t BjB k ) (1 - f ijk ) $ 2 (^ - ipj;v m (i), v m (j)) + (jki) 

- (2^)" 3 (BiBjB k ) f ijk + (jki) 

+ constants (22) 

This cumulant includes components typical of jets (line 1 and 2), flow (line 4-11'), and one term unique to jet-flow 
cross-term (line 3). While this jet-flow term does include two-particle jet-like, and flow-like factors, it is important to 
note the flow factor is intrinsically inharmonic. Its amplitude depends on the number jet associates and background 
particle multiplicity. Given the number of associates is likely much smaller than the number of background particles, 
one expects the amplitude of this term to be dominated by the magnitude of the background. It is interesting to 
compare the magnitude of this cross-term relative to the off-diagonal flow terms. We focus on the difference between 
the leading terms in v 2 v 2 V4 and v 2 v 2 . Given the jet cross section and fragment multiplicity are small at RHIC energies, 
one expects the v 2 v 2 V4 term should dominate or at most be of similar magnitude to the cross-term unless the jet flow 
v 2 is considerably larger than the bulk flow. 

As discussed in Sect IIII 51 particle emission in A + A collisions is a non-Poissonian process. The coefficients fij k in 
general deviate from unity. The leading background terms may thus be those of line 2, 6, 7, and 9 in Eq. 1221 

Lines 4 and 7 both contain terms in F 2 albeit with different amplitudes. Experimentally these cannot be distin- 
guished, and one ends up with an F 2 contribution to the cumulant which depends on both the jet and background 
yields. Likewise, lines 5, 6, 8, and 9 feature non-diagonal flow terms (dominated by v 2 v 2 v^) which cannot be distin- 
guished experimentally: the amplitude of the v 2 v 2 v^ three particle correlation terms depend intricately on the jet yield, 
its fluctuations, as well as the background yield. The jet-flow 'cross-term' manifests itself through the inharmonic 
terms of line 3 in Eq. [22j While these harmonic terms may be approximated as a product such as G2(A</>i2)F<I>2(A</>i3) 
if the width <j\ (very high-p t particles) is negligible, when added in quadrature, to the width a 2 of the low-p t particle, 
the approximation breaks down in general because o\ is neither zero or equal to a 2 . Additionally, note that the model 
used in this section assumed the cross-terms are dominated by jetty and flow components. This may not be the case 
in practice. For instance resonance or cluster decays in the presence of both radial and elliptical flow should lead to 
complex cross terms. Such cross-terms shall be inhenrently inharmonic also. Given the correlation shapes associated 
with resonance or cluster decays can be quite wide, the inharmonic character of the cross term can be quite intricate. 
Adhoc subtraction of terms in G2(A</>i2)F<i>2(A</>i3) in model based analyses is therefore unwarranted. 

IV. PROBABILITY CUMULANTS 

We showed in Section IIII Bl that the 3-cumulant associated with flow reduces to non- diagonal terms in v n v m v p 
for Poissonian particle production processes. However, momentum, energy, and quantum number conservation laws 
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imply elementary collisions are typically non-Poissonian processes as observed from Fig 3 which shows ratios fijk 
measured by STAR (46j . This implies the 3-cumulant may have a rather complicated structure, with second order 
terms (i.e. two-particles) as well as three particle terms. We note however that the simplicity of the 3-cumulant may 
be recovered, in this case, by using probability cumulants rather than density cumulants. The probability cumulants 
are defined as follows: 

P2( ^ } = - m (n s ) (23) 



P3(<Pi,<Pj,<Pk) P2((fii,<Pj)pi{<Pk) 



(NiNjNk) {NiNj) (N k ) 

(NiN k ) (Nj) (NjN k ) {Nt} 

2 Pi(Vi)Pi( l Pj)Pi( l Pk) 
(N t ) (Nj) (N k ) 



(24) 



where JVi correspond to total particle multiplicities accepted in the kinematic cuts i. It is straightforward to verify 
that the flow probability 3-cumulant indeed consists of non-diagonal flow terms only: 

Pi , (cp 1 ,(p 2 ,( P3 ) = (2Tr)- s ^ 3 ((p 1 ,( P 2,cp3;v p (l)v m (2)v n (3)) (25) 

The probability cumulant, P3, thus provides a tool to suppress the strength of non-Poissonian 2 nd order terms, and 
may therefore be used, in addition to the number 3-cumulant, C3, in three-particle analyses. 

The simplification obtained for flow processes is unfortunately not realized for jet or conical emission processes. 
The use of probability 3-cumulant may then be of limited interest in practice. This can be straightforwardly shown 
through a calculation of the probability cumulant of the Gaussian model used in Sec. IIII A[ We limit our calculation 
to the near side jet. One finds the probability 3-cumulant contains terms in G 2 (A(f>) proportional to the non vanishing 
factor (J(J - 1)) (NiN 2 ) (N 3 ) - (J) 2 (NiN 2 N 3 ), where Ni, N X N 2 , and N ± N 2 N 3 are respectively the number of singles, 
pairs, and triplets from the near side jet. By contrast to the flow probability 3-cumulant which contains only genuine 
three-particle correlation terms ($3), the Gaussian probability cumulant thus has a complicated expression which 
contains terms in G3 as well as in G 2 . Its interpretation is therefore non trivial. 

V. EFFICIENCY CORRECTION AND OBSERVABLE ROBUSTNESS 

Measured particle densities, and cumulants must be corrected for detection efficiencies and other instrumental 
effects. We introduced in Ref. [4l[ a procedure to correct for detector efficiencies based on ratios of two- and three- 
particle densities by products of two, and three single particle densities. For instance, up to a global efficiency factor, 
the corrected three-particle density is given by 

P3(<P1, <P2, ^corrected = N(1)N(2)N(3) N * ^ ' f\ (26) 

^1(^1)^1(^2)^1(^3) 

where Ni(tpi) stands for uncorrected single particle densities, N 3 (ipi, ip 2l f 3 ) are numbers of triplets, and N(i) are 
total particle yields within the kinematics cuts i. 

This correction procedure is strictly exact for continuous functions provided the efficiency is a function of the 
azimuthal angles (fi but independent of other observables such as the particle rapidity and transverse momentum. For 
large detectors such as the STAR TPC |4j|, the efficiency is a smooth and slowly varying function of the pseudorapidity 
and transverse momentum, but exhibit periodic structures in azimuth because of TPC sector boundaries. A correction 
for azimuthal dependencies of the detector response is thus the most important and relevant in this context. 

We examine the robustness of the above correction procedure in a practical situation, i.e. where the densities 
are measured with a finite number of bins, based on the jet and conical emission models discussed in Sec. 2. We 
parameterize the detection efficiency, £(<fi), with a Fourier decomposition. 

e(tpi) =e ^1 + 2 Y^ £ i cos ( ni Pi)\ (27) 

where eo is the average efficiency, and coefficients e„ determine the azimuthal dependence of the efficiency. We 
assume the efficiency for simultaneously detecting two and three particles arc factorizablc as product of single particle 
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efficiencies. The mean single particle detection efficiency is set to 80%. The Fourier coefficients (Eq. [27]) are chosen 
to obtain a non-uniform azimuthal detector response as shown in Fig. [3J The JET+Mach Cone model described in 




cp (rad) 

FIG. 3: (Color online) Detection efficiency dependence on the particle azimuthal angle used in the simulation results presented 
in Fig. H 

sections IlII Cl and lVIl is used to carry a simulation of the robustness of the correction procedure. Figure |4] (a) displays 
the 3-cumulants obtain with perfect detection efficiency. Figure U (b) shows the uncorrected two particle density 
Pii^x^n) obtained with the non-uniform azimuthal response illustrated in Fig[3l The two particle density exhibits 
strong, narrow, and repeated structures that result from the twelve-fold structure of the assumed detector response. 
Figure 0] (c) displays the 3-cumulant obtained with the same detector response, but corrected for efficiency effects 
by division of the two and three particle densities by the product of single particle densities, as in Eq. 1261 In stark 
contrast to the correlation shown in Figure[4](b), one finds the corrected cumulant exhibits no evidence of the detector 
response except for a global change in amplitude corresponding to the cube of the average detection efficiency. Figure 
S] (d) shows the difference between this cumulant and the perfect efficiency cumulant scaled down by the cube of the 
average detection efficiency. The standard deviation of the difference is smaller than 1% of the maximum 3-cumulant 
amplitude and exhibit no particular structure. We thus conclude the correction method is numerically robust for 
cases where the triplet, and pair efficiencies factorize. 

VI. CUMULANT METHOD SENSITIVITY FOR JET AND CONICAL EMISSION MEASUREMENTS 

We present in this section a study of the sensitivity of the cumulant method described in this paper for a search 
for conical emission. Our study of the sensitivity of the method is based on Monte Carlo (MC) simulations carried 
out using a simple event generator that encapsulates flow, jet, and conical emission as a multicomponent model. 

We begin our presentation of the simulations with flow-only model components. The overall multiplicity, noted 
m, of an event is selected randomly according to a flat distribution. The flow component is designed to produce, 
on average, a fraction, fi, of the event multiplicity. The number of flowing particles, iVj, (for a kinematic cut i) is 
generated event-by-event randomly with a Poisson probability density function (PDF) of mean m x The particle 
direction, in azimuth, is selected randomly according to the flow PDF given by Eq. [7] Various values of the coefficients 
i>2 and V4, which determine the magnitude of the bulk flow, are used in the following to study the impact of flow on 
the cumulant and sensitivity to a cone signal. All other Fourier coefficients are set to zero. The event plane azimuthal 
angle is chosen randomly in the range [0, 2tt[. The coefficient /i is set to produce a number of high-p f particles or order 
unity in each event, while /2 is set to produce a low p t particle multiplicity of order 100. These values are selected to 
mimic STAR data in the range 3 < pt < 20, and 1 < pt < 2 for the trigger and associate particles respectively. 

Figure[5](a) and (b) respectively display the triplet density and the 3-cumulant obtained with the above flow random 
generator with v-2 = 0.1, and V4 = 0. The 3-cumulant, shown in Fig[5](b), exhibits a finite v\ component owing to the 
fact that although the multiplicities Ni are generated with Poisson PDFs, the mean of these PDFs varies as a function 
of the event multiplicity thereby implying a correlation between the low and high p t particle multiplicities. This, in 
turn, implies the coefficients fijk are non-zero: the cumulant therefore contains a non-vanishing v\ component. This 
component however vanishes in the probability cumulant shown in Fig [5] (c) as expected from Eq. 1251 

Figure [5] shows results of a simulation based on the same flow model as that used in Fig. [3J but with an added 
fourth harmonic component. Flow amplitudes are set to V2 — 0.1 and V4 = 0.01. While the 3-particle density (Fig [(5] 
(a)) is completely dominated by the 2 nd harmonic, one finds the 3-cumulant suppresses this component significantly 
(Fig [5] (b)), and enables clear observation of the v\v± non-diagonal terms expected from Eq. [2U The v\ component 
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FIG. 4: Study of the efficiency correction method (Eq. I26(l using the Mach Cone model described in Sec. IVII (a) 3-cumulant 
for perfect detection efficiency, (b) Two particle density, (c) 3-cumulant obtained with the azimuthal efficiency dependence 
shown in Fig. [3] (d) Difference between the 3-cumulants shown in (b), and the perfect efficiency cumulant scaled by the cube 
of the average detection efficiency. 



is further suppressed in the probability cumulant, shown in Fig [5] (c), where only the v\v^ component manifestly 
remains. As anticipated based on Eq. I25i the probability cumulant enables essentially full suppression of v\ terms, 
and show irreducible flow components only. 

Mono-jets and conical emission are next added to the simulated events. The jet axis (corresponding to the parton 
direction initiating the jet) is chosen randomly in the transverse plane. Particles are generated at random polar angles, 
relative to the jet direction using Gaussian PDFs. The width of the Gaussian is set to 0.15 radian for high p t particles, 
and 0.25 for low p t particles. Note that the conclusions of this study are essentially independent of the width of the 
jets. The associate multiplicity are generated jet-by-jet using Poisson PDFs. The near side jet associated particle 
multiplicity (i.e. number of associates per jet) are set to 1 and 2 for the high- and low-p t particles respectively. No 
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FIG. 5: Simulations of (a) the three-particle density, (b) three-cumulant, and (c) the probability cumulant for particle emission 
with elliptic anisotropy V2 — 0.10. 



(a) 



to 



9- 
< 





00 




Acp i2 (rad) 



A(p i2 (rad) 



A(p i2 (rad) 



FIG. 6: Simulations of (a) the three-particle density, (b) three-cumulant, and (c) the probability cumulant for particle emission 
with elliptic anisotropy V2 — 0.10, and V4 = 0.01. 



away side jet is included. Instead, one introduces a conical emission as described in the following. Jet parameters 
are selected to correspond approximately to associate jet yields measured by the CDF collaboration for jets of 10 to 
20 GeV [13], as well as jet associated yields reported by RHIC experiments [HI]. The number of jets is generated 
event-by-event with a Poisson PDF of mean 0.01 jet per unit multiplicity: the average number of jets per event is 
thence of order 1-2. Such large values are used to ease our study. Smaller, perhaps more realistic jet-yield values result 
in a weaker signal. This limits the statistical significance of the measurement, but does not constitute an intrinsic 
limitation of the cumulant method. 

In our simulations, the parton determining the direction of the jet is first generated, and is restricted for simplicity 
to the reaction transverse plane. The away-side parton is emitted in the transverse plane as well, at angle of tt radian 
relative to the jet parton. Conical emission is simulated by generating particles at a random polar angle relative to the 
away-side parton direction. The azimuthal direction of the particle, in the plane transverse to the parton direction, is 
randomly selected based on a uniform distribution. The polar angle is chosen to have a Gaussian distribution with an 
average value of one (1) radian. The number of high-p t particles produced in this cone is arbitrarily set to zero, while 
the number of low-pt cone particles in a given parton-parton collision is selected randomly with a Poisson distribution 
of mean equal two, unless specified otherwise in the following. 

Because the cone particles are emitted at an arbitrary polar angle, but detected in the transverse plane as a 
function of relative azimuthal angles only, conical emission results in four Jacobian peaks, as seen in Fig. [7| (a). In 
order to characterize the shape and strength of these peaks, we use projections of the cumulants along the diagonals 
V12 + ¥>13 — 2-7T and ipi2 — fn- Projections, shown in Fig. [7] (b) are limited to include narrow angular ranges about 
the diagonals. Specifically, the range \tp\% — y>\z\ < 10° is integrated to obtain the projection along the plot main 
diagonal (displayed in red), whereas the condition \((pi2 + V13) — 27r| < 10° is used for projections along (p±2 — <pia 
(shown in blue). 
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At issue is whether or not one can extract a cone signal given the presence of background particles, and flow. We 
thus vary the strength of the jet signal to background ratio by changing the fractional jet and background yields. 
Figures [7] through [Tl] show cumulants and projections obtained with various model parameters. All plotted cumulants 
were obtained with simulations integrating two millions events. 

Figures[jj(a-b) shows the cumulant and projections obtained for an average of 2.5 jets per event. The number of high 
and low pt particles associated with the jet, and cone are set respectively to 1 and 2. The uncorrelated background 
of high- and low-p f particle are set to 2.5 and 100 particles per event respectively. The cone signal obtained with 
such parameters is obviously strong and clearly observable. Figure [7] (c) displays the cumulant and projections with 
a high-p t background raised by a factor of 4 resuling in a signal to noise ratio of 25%. The observed cone strength 
remains unchanged, as expected given the background is uncorrelated. We also verified that an increase of the width 
of the gaussian PDF used for jet, and cone particle emission results in wider peak structures with reduced amplitude 
but no actual change in the integrated cone signal strenght. In this context, we conclude the observability of conical 
emission is only limited by the statistical accuracy of the measurement relative to the actual strength of the signal. 
For instantce, a reduction by a factor of five of the number of high-p t associated particles results in cone-signal five 
times small in amplitude, but the shape of the cumulant remains unchanged and the observability of the signal is thus 
only limited by the size of the event sample, and the number of particles associated with the trigger and cone. 
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FIG. 7: (a) Simulation of the 3-cumulant obtained with mono-jet, conical emission, and background particles (without flow). 
Average jet yield to high-pt particle background yield set to 1. (b) Projections of the cumulant along the main (red) and 
alternate (blue) diagonals for jet/cone parameters used in (a). Other particle generation parameters, and projections are 
discussed in the text, (c) Cumulant projections obtained for a high-pt background to jet-signal ratio of 4. 



We next consider the observability of conical emission in the presence of anisotropic flow. Figure [8] shows the 
cumulant and projections in the presence of a flow background. The jet and conical signals are identical to those used 
in Fig. [7J (a-b). An average background of 100 low-p t particles is included with v-2 = 0.1, and V4 = uf. No high-p t 
background is included. The cone signal shape and strength remains unaffected, although the projections exhibit 
a background structure somewhat more complex than that observed in Fig. [Jj Fig. [8] (c) displays the projections 
obtained when the elliptic flow is raised to i>2 = 0.2. While the signal is clearly visible, one observes the emergences of 
cosine structures along the Ay?23 projection. These stem from the non-Poisson nature of the background multiplicity 
fluctuations used in the simulations. 

Flow-like or cosine structures also appear when high-p t background particles are added as shown in Figure [5] where 
the average high-p t particle background is set to 2.5/event. The observability of the cone depends on its amplitude 
relative to that of the background, as well as on the magnitude of V2 and V4. We note that the signal remains visible 
in the Ac^23 projection, with the jet and cone parameters used in Figure [JJ even when the high-p t background is 
raised by a factor of eight as shown in Fig. [9] For substantially larger background, the signal however becomes more 
difficult to extract as illustrated in Figure [TO] where a background to signal ratio of sixteen is used for high-p t particles. 
Clearly, the observability of a cone signal in the presence of flow depends on the strength of the signal relative to that 
of non-Poisson components and values of i>2- 

The model independent extraction of a conical signal becomes increasingly difficult as the signal strength is reduced 
and the elliptic flow magnitude increased. The measurement becomes particularly challenging when the jet amplitude 
is modulated by differential attenuation relative to the reaction plane, i.e. in the presence of jet-flow. Figures [Tl] (a) 
and (b) display projections along Aipi2 + A(pi3 and A</?23 obtained with jet V2 — 0.05, (as defined in Sec. IIII D[) . 
The blue squares and red circles show, respectively, the signals obtained with an average of two and one low-p t 
particles involved in conical emission. Simulations are carried out with an average number of 2.5 jets/cones per event. 
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FIG. 8: Idem Fig. [7] in the presence of background flow, (a-b) No background high-p t particle. Average of 100 low-p t 
background particle per event, with flow V2 = 0.1, and v± = wf. (c) Average of 100 low-p t background particle per event, with 
flow V2 — 0.2, and V4 = v\. 
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FIG. 9: (a-b) Idem Fig. [7] for average jet yield to high-p t background equal one per event. Average of 100 low-pt background 
particle per event, with flow V2 = 0.1, and v± = v\. (c) Projections obtained for average jet yield to high-p t background equal 
eight (8) per event. 
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FIG. 10: Idem Fig. [Jjfor average high-pt background to jet yield equal sixteen (16). Average of 100 low-pt background particle 
per event, with flow V2 = 0.1, and V4 — v\. 



The average yield of high-p t low-pt background particle are set to 8 and 100 per event respectively. Background 
particles are produced with flow vi = 0.1, and v± = uf. In our simulation, the number of cone associated particles is 
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determined event-by-event with a Poisson random number generator. The conical signal is detectable, in the context 
of this 3-cumulant analysis, only when two, or more, particles are generated per event. For an average cone yield of 
two associated particles, 60% of the cone events contain two or more cone particles. By contrast, only 26% of the 
events contain two cone particle when the average (cone) yield is one. The signal obtained in the former case is clearly 
visible, while if becomes weaker in the latter. The detectability of the signal thus greatly depends of the actual yield 
of conical emission. We note that two-particle correlation analyses report an away side yield of two particle per trigger 
after flow background subtraction based on a ZYAM hypothesis [5l[ . If this away side is due to conical emission, this 
implies the \ow-p t cone particle average multiplicity is large and should therefore be detectable in the context of a 
3-cumulant analysis. We note however that the cumulant measurements reported by the STAR collaboration at recent 
conferences do not exhibit as strong cos(2A(p) components as those illustrated in the Atpi + Aip 2 projections shown 
in Fig. [11] [42||. Note additionally that the cone signal remains visible in the Aip 2 3 projection even for small signal 
to background ratios. It is however obvious that strong background flow and jet-flow may hinder the observation 
of conical emission. For instance, Figure [12] shows a comparison of the 3-cumulant projections obtained for conical 
emission in the presence of jet-flow cross terms with V2 values of 0.05 (blue) and 0.1 (red) an associate cone yield of 
one. One finds the cone single is difficult to distinguish for vi = 0.1 but straightforward to observe for V2 = 0.05. 

In summary, e find in the context of our multicomponcnt model that the exact sensitivity of the measurement is 
a " complicated" function of the flow and background parameters of the model, and cannot readily be expressed in 
analytical form or even as a table. We further note that the shape and strength of actual correlations is influenced 
by momentum, and quantum number conservation effects not explicitly addressed in this work. 
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FIG. 11: (Color online) Comparison of the cone signal obtained with an average of two (blue squares) and one (red circles) 
low-pi particles involved in conical emission. Panel (a) and (b) show the 3-cumulant projections along (A(p%2 + Ay>i3)/2 
and Aif23 respectively. Simulations carried out with a jet-flow component, V2 = 0.05, and an average number of 100 low-p t 
background particles per event, with flow v% = 0.1, and V4 = v\. See text for details. 



VII. SUMMARY AND CONCLUSIONS 



Conical emission has been suggested in the recent literature to explain the away-side structures observed in two 
particle azimuthal correlations for Au + Au collisions at ^/s n n = 200 GeV. If realized, conical emission should also 
lead to three particle correlations. We introduced in Ref. [41j a method based on cumulants to carry a search for 
such a signal. In this work, we presented a discussion of the properties of the three particle cumulant (3-cumulant) 
calculated in terms of two relative azimuthal angles, as defined in [4l|, for searches of conical emission in high-energy 
A+ A collisions. We showed that in the presence of azimuthal anisotry (flow), the 3- cumulant reduces to non-diagonal 
terms dominated by components of order V2V2V4 for Poissonian particle production. Given particle production is in 
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FIG. 12: (Color online) Comparison of the cone signal obtained with an average of one low-pt particles involved in conical 
emission and an average of 100 low-pt background particles per event with V2 = 0.05 (red circles) and V2 = 0.1 (blue squares). 
Panel (a) and (b) show the 3-cumulant projections along (A<^i2 + A<^i3)/2 and A<^23 respectively. Simulations are carried out 
with a jet-flow component, V2 = 0.05. The background has V4 = v\. 



general non-Poissonian, one however expects the presence of second order terms, dominated by components in v\ in 
the cumulant. The strength of those terms relative to the irreducible V2V2V4 components depends on the magnitude 
of V2 and the global strength of particle correlations, or departure from Poisson statistics. The strength of particle 
correlation is known to vary inversely with the collision multiplicity in A + A collisions due to increasing two-particle 
correlation dilution with increasing number of collision participants. We thus expect the mix of V2V2V4 and v\ should 
be a function of collision centrality. 

We introduced a three-particle probability cumulant, and showed it is devoid of v\ non-Poissonian components in 
the presence of flow only thereby enabling, in principle, the determination of V2V2V4 amplitudes. We discussed the 
shape and strength of the 3-cumulant based on simple particle production models including di-jets, mono-jet plus 
conical emission, and jet x flow correlations. We showed that jet x flow cross-terms should be more complex than those 
assumed by some ongoing analyses, and cannot, in particular, be derived from a product of two-particle flow and 
jet-like terms. We used the models to discuss the sensitivity of the 3-cumulant for searches of conical emission in 
A + A collisions. We showed that the cumulant enables excellent sensitivity to conical emission for modest values of 
flow i>2, but becomes increasingly challenging for larger values of V2- We further showed that conical emission can be 
detected even in the presence of jet x flow correlations for small values of i>2 jet flow. It is however obvious the signal 
may be masked by large values of jet and background V2 and large non-Poissonian particle production. Precise values 
of sensitivity, in terms of signal strength to background ratios, cannot be expressed in a model independent way, and 
are found to depend on the relative amplitude of background V2 and i>4 jet-flow correlations, as well as the global 
strength of particle correlations, i.e. departure from Poisson statistics. 

This work neglects effects associated with momentum conservation, quantum number conservation, and particle 
correlations present inp+p collisions. Momentum conservation implies in particular that a multicomponent description 
of particle production is not strictly valid, and while it is useful to estimate the effects of various particle production 
mechanisms, its use to search for conical emission signals is model dependent and therefore unreliable. 
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